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Abstract 

Recently we have proposed a simple hadronic axion model within gauge- 
mediated supersymmetry breaking. In this paper we discuss various cosmo- 
logical consequences of the model in great detail. A particular attention is 
paid to a saxion, a scalar partner of an axion, which is produced as a coherent 
oscillation in the early universe. We show that our model is cosmologically 
viable, if the reheating temperature of inflation is sufficiently low. We also 
discuss the late decay of the saxion which gives a preferable power spectrum 
of the density fluctuation in the standard cold dark matter model when com- 
pared with the observation. 
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I. INTRODUCTION 



The most attractive candidate for the solution of the strong CP problem is the Peccei- 
Quinn (PQ) mechanism [|lj. In this mechanism, there exists an axion which is the Nambu- 
Goldstone (NG) boson associated with the global U(1)pq PQ symmetry breaking. 

In the framework of gauge mediated supersymmetry (SUSY) breaking theories M, we 
have proposed an interesting possibility || to dynamically generate the PQ symmetry break- 
ing scale in a so-called hadronic axion model M. A gauge singlet PQ multiplet X and colored 
PQ quark multiplets Qp and Q P are introduced with the superpotential 

W = \pXQ P Qp~, (1) 

where Xp is a coupling constant. Here we stress that no mass parameter is introduced. Be- 
cause in a supersymmetric limit the U(1)pq symmetry is enhanced to its complex extension, 
there appears a flat direction Q P = Q P = with X undermined in the same limit. The 
effects of the SUSY breaking, however, lift the flat direction and fix the vacuum expectation 
value (vev) of the scalar PQ field X. In Ref. 0] we have shown that the balance of the 
SUSY breaking effects between the gravity mediation and the gauge mediation stabilizes X 
and gives the nonzero vev < X >= fpq approximately as 




where m a denotes the mass of a saxion field, the real part of the scalar component of X. We 
estimate m s = ^3/2 with a parameter £ of order unity. And / is a mass scale fixed by the 
messenger sector of the gauge meditation mechanising A current scalar lepton mass limit 
suggests that / > 10 4 GeV. In the following analysis, we take the minimum value / = 10 4 
GeV for a representative value. 

1 / is roughly represented as f 2 ~ (f*) {F), where (F) denotes the F-component vev of the 
messenger singlet field and a s is the QCD coupling. 
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The axion arises when the X field develops the non-vanishing vev. The decay constant 
of the axion (the PQ scale) fpQ is constrained by various astrophysical and cosmo logical 
considerations. In particular the cooling of the SN 1987A puts a lower bound as fpQ ^ 10 9 
GeV [||]. On the other hand, the upper bound on Jpq comes from the overclosure limit of 
the axion abundance. The relic abundance of the axion due to the vacuum misalignment is 
given by f||7| 

with 9 being an initial misalignment angle (\6\ < ir) and h the present Hubble constant in 
units of 100 km/sec/Mpc. Thus this leads to fpQ ^5 10 12 GeV. Here it has been assumed 
that no entropy is produced after the QCD phase transition. If extra entropy production 
mechanism with a reheat temperature Tp is operated after the transition, the estimate of 
the relic abundance changes as |||9j 

^ 2 -^ 2 (lotv)(lo^v) 2 - <«> 

Then the upper bound on the PQ scale becomes fpq < 10 16 GeV for T R ~ 10 MeV. Note 
that the big bang nucleosynthesis implies that the reheating temperature should be higher 
than about 10 MeV. 

It follows from Eq. (Q) that, in the model we are considering, the allowed region for the 
PQ scale discussed above 

10 9 GeV < f PQ < 10 12 (10 16 ) GeV (5) 

can be translated to the allowed region for the saxion mass 

100 MeV > m a ~ m 3/2 > 100(10~ 2 ) keV. (6) 

Remarkably many models of the GMSB predict a gravitino mass in the above range. Thus 
the model gives a very simple description of the PQ breaking mechanism solely governed by 
the physics of the SUSY breaking. 
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In the model, however, the saxion has a very light mass which is comparable to the 
gravitino mass and its lifetime may be long. Furthermore the saxion might be produced as 
a coherent oscillation too much. Therefore, in this paper we will consider the cosmology of 
the model, especially the saxion cosmology, in detail. 

The paper is organized as follows. In the subsequent section, we will discuss cosmological 
evolution of the saxion and estimate its relic abundance. In Section II, we will consider 
cosmological constraints on the saxion and show that our model survives the constraints if 
the reheating temperature of the inflation is sufficiently low. Then in Section IV, we will 
point out the possibility that the saxion plays a role of a late decaying particle, giving a 
better fit of the density perturbation in the cold dark matter scenario. Final section is 
devoted to conclusions. 



II. COSMOLOGICAL EVOLUTION OF SAXION 

Let us begin by discussing the cosmological evolution of the saxion particle in an inflation- 
ary universe, and estimate its relic abundance. Since the saxion has interaction suppressed 
by the energy scale Jpq, it decouples from the thermal bath of the universe when the cosmic 
temperature becomes about pi 



r - ~ 108 GeV lira) ~ 109 GeV (re) (^r) • < 7 > 

If the temperature achieved after the primordial inflation was higher than T^ ec , the saxion 
could be thermalized and the ratio of the energy density of the saxion p a to the entropy 
density s is estimated as 



10" 6 GeV — — . (8) 
s VI MeV/ V ; 

Even if the saxion was not thermalized, the saxion could be produced by scattering processes 
right after the reheating era. In this case the saxion abundance {pa/s) is given by 

^ ~ 10- 6 GeV ( m ° ^( Tr 



s VI MeV7 \T dec 



-lO^GeV f-gyn / 1 f Tfl ), (9) 
VIMeVy V10 4 GeV; UO MeV 7 ' V; 

where Tr is the reheating temperature of the inflation (Tr > 10 MeV) and it is highly 
model dependent [|11|]. In the gauge-mediated SUSY breaking models, since the longitudinal 



component of the light gravitino has interaction much stronger than the gravitational one, 
it would be abundant and overdose the universe unless the mass is lower than about 1 keV. 
To cure this, the reheating temperature Tr must be sufficiently low |T^JT3|] p| 



m~ /VI MeV 

With such a low Tr, the relic abundance of the saxion produced by the scattering in the 
thermal bath is highly suppressed. 

On the other hand, the saxion is also produced as a coherent oscillation. The scalar 
potential in general receives corrections at an early universe, e.g. during the primordial 
inflation. Thus the saxion does not sit at the true minimum a = Jpq but is generically 
displaced at o = 07 after the inflation. Knowledge of the initial displacement 07 is crucial to 
estimate the abundance of the oscillating saxion, but it depends on the details of supergravity 
Lagrangian as well as inflation models.^ Thus in the following, we regard 07 as an arbitrary 
parameter which is not larger than the Planck scale (Mq = 2.4 x 10 18 GeV) and consider 
two distinct cases. 



2 If one assumes that the universe experienced the late-time entropy production with a sufficiently 
low reheating temperature, the gravitino is sufficiently diluted and we have no upper bound on the 
reheating temperature of the primordial inflation Tr like Eq.(10). 

3 07 ~ or 07 ~ Mq may be naturally expected by the supergravity effect. 
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A. Case 1: Saxion trapped at the origin 



Let us first consider the case that the saxion field sits around the origin (07 ~ 0) just 
after the inflation.^ In this case the PQ quark multiplets (Qp and Qp) are in the thermal 
equilibrium through the rapid QCD interaction, which generates an additional mass term 
~ T 2 |Y| 2 to the saxion potential as a finite temperature effect. Thus the saxion is trapped at 
the origin till the cosmic temperature becomes below a critical temperature Tc- We expect 
that Tc ~ 100 GeV, since the saxion obtains a negative soft SUSY breaking mass squared 
at the origin due to the renormalization group effect of the Yukawa interaction Eq.([IJ) and 
its absolute value is expected to be of the order of the electroweak scale ~ 100 GeV. 

Therefore the saxion starts to roll down toward its true minimum at T = Tc and oscillates 
around it with the initial amplitude cr ~ /pq-| The ratio of the energy density of this 
oscillation to the entropy density is estimated as 

1 2 f 2 ( a o y 

where g* counts the degrees of freedom of the relativistic particles (g* ~ 100 at T = Tc). 
Here it should be noted that this ratio takes a constant value till the saxion decays if no 
extra entropy is produced, since p a and s are both diluted by p a oc s oc R~ 3 (R: the scale 



4 Thus the U{1)pq symmetry is recovered. In this case it should be noted that, since the effective 
mass which the saxion feels during the primordial inflation is generally comparable to the Hubble 
parameter during the inflation, the quantum fluctuations for the X field are strongly suppressed 



[14]. For special Kahler potential, a massless mode appears in the X field and a large quantum 
fluctuation is expected. This would lead to too much axion isocurvature fluctuation. To suppress 
it, the mass density of the axion or the Hubble parameter during the inflation should be small. We 
thank M. Kawasaki and T. Yanagida for pointing out this problem. 
5 We are free from the domain wall problem in our model, since Np>y/ = 1. 
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factor of the universe) as the universe expands. In the model we are discussing, the PQ 
scale fpQ is given by Eq.(^j) and thus the abundance of the saxion oscillation reads 



^ ~ 1.1 x 1(T 2 ^ ( a ° 



2 



s ' n \f PQ ) ' 

- llxlo8GeV (lo4v) 4 (lo^v)"fe) 2 ' < 12 > 

Note that it is independent of the saxion mass. Here we have assumed that the reheating 
process of the primordial inflation completed before the saxion oscillation starts, i.e., Tr > 
T c . On the other hand, for the case Tr < T c , we expect a dilution of the saxion energy 
density. The saxion energy density at T = Tr is given by 

^-^)(li) 3 ^^fe) 2 @i) 3 - (13) 

From the fact that the radiation energy density is proportional to R~ 3 ^ 2 from the end of the 
inflation till T = Tr, it follows that the ratio of the scale factor (R(Tc) / R(Tr)) becomes 

(R(T C )\ 3/2 g*(T R )T R 



\R(T R )J g.(T c )T%' 
Thus the saxion abundance at T = Tr is given by 



(14) 



1.2xl(T 13 GeV ( —J— V (-J°—)~*(-J*^- V ( V . (15) 
[ 10 A GeV UOOGeW UO MeVj 1 ,: ' 



Here note that the reheating temperature Tr should be higher than about 10 MeV in order 
to maintain the success of the big bang nucleosynthesis. Comparing it with Eq. (|I~2"D , we find 
that the relic abundance is diluted significantly by the primordial inflation with a very low 
Tr. 

The above discussion does also hold even for the case where the saxion is displaced from 
the origin 07 by a small amount. This is because the PQ multiplets will be thermalized and 
generate the potential which traps the saxion at the origin, as far as their effective mass 
(~ ApOf) is smaller than the maximum temperature achieved after the primordial inflation. 



6 The maximum temperature Tmax achieved after the primordial inflation is given by |?J as 
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B. Case 2: Saxion displaced far from the origin 



If the displacement 07 is large enough, the cosmological evolution of the saxion is com- 
pletely different. Here, we consider only the case that (M G ^)oj > /pq-0 When the Hubble 
parameter of the universe (H) becomes H ~ m CT , the saxion starts to obey the coherent 



Using the initial amplitude of the oscillation cr = 07 (> /rq), the saxion abundance of this 
oscillation is estimated as 



Here we have assumed that the saxion oscillation starts after the reheating of the primordial 
inflation had been completed (Tr > T a ). However, such a high Tr leads the gravitino 
problem for m a (~ 7713/2) ^ 1 keV. Comparison of Eq. (|i~6|) with the upper bound on Tr 
[Eq.fllOD] shows that T CT is always higher than this bound. On the other hand, the inflation 
with Tr < T a dilutes the energy density of the saxion as 



Tmax 0.31% Vj ' , where V 1 is the vacuum energy of the inflation. 

7 Thus the U{\)pq symmetry is broken and a massless mode, i.e. the axion, appears. See the 
footnote 4 for an argument on the density fluctuation of the axion field. 

8 Note that for the region a = aj /pq the potential is governed by the contribution from the 
gravity mediation effect, and the effective mass of the saxion is almost the same as the gravitino 
mass, i.e., the saxion mass (m 3 / 2 — m^). 



oscillation.[| The temperature at this time T a is estimated as 




(16) 




(17) 



S 



s ~ 8 R \M, 



G 

~ 1.3 x 10- 3 GeV ( (-£-)*. (18) 
V10 MeV/ \MqJ v ; 

This ratio takes its minimum value for the lowest reheating temperature Tr ~ 10 MeV and 
the smallest initial amplitude cr ~ f P g. 

To end this section we estimate the saxion lifetime. When the Hubble parameter becomes 
comparable to its decay width, the saxion decay occurs. In the present model, the saxion 
dominantly decays into two axions with the ratef] 

(19) 



and the lifetime is given as 

^"""-WdMr- <2o) 

Therefore the saxion becomes stable within the age of the universe (~ 10 17 sec) if its mass is 
smaller than about 10 keV. Note that the value r a varies very widely because of the strong 
dependence on the saxion mass. 



III. CONSTRAINTS ON COSMIC SAXION 

In this section we consider the cosmological constraints on the saxion abundance and 
examine whether the proposed model is viable or not. Since the energy density of the saxion 
is not so suppressed compared to the radiation and also it may have a long lifetime, we have 
to carefully see its cosmological effects. 

As described in the previous section, the saxion is produced by the scattering at the 
reheating era and also as the coherent oscillation. Comparing the both saxion abundances 



9 The saxion may decay into two photons with a very small branching ratio. We will discuss this 
rare process in the next section. 
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with each other, one finds that the saxion produced from the thermal bath can be neglected, 
if one assumes the primordial inflation with a relatively low reheating temperature to avoid 
the gravitino problem and too much relic abundance of the saxion oscillation. Therefore it 
should be enough to consider only the saxion abundance coming from the oscillation in the 
subsequent discussion. 

A. Constraints on Saxion Abundance 

First of all, the saxion with r CT > 1 sec, i.e., m a < 100 MeV, may upset the big bang 
nucleosynthesis. Existence of the extra energy of exotic particle at T ~ 1 MeV may speed up 
the expansion of the universe and increase the number ratio of neutron to proton, resulting 
in the overproduction of 4 He. Roughly speaking, the energy of such a particle should be 
smaller than that of one neutrino species. Thus the saxion abundance should satisfy 

P — < 10- 4 GeV. (21) 

s 

Furthermore, if the saxion lifetime is larger than the age of the universe (i.e., m a < 
10 keV), the saxion oscillation still exists now. Not to overdose the universe, the saxion 
abundance should be 

El < El = 3.6 x 10- 9 h 2 GeV, (22) 
s s 

with the critical density p c and the present entropy density s . 

Next, we turn to the rare decay of the saxion. Since the saxion is light enough in our 
model, it can not decay into two gluons but decays into two photons via one-loop diagrams 
if the PQ quark multiplets have QED charges. Thus we expect that the branching ratio of 
this radiative decay is very suppressed as B 1 ~ (a em / (47r)) 2 ~ 10~ 7 . For example, if the PQ 
quark fields Qp and Q P are in 5 and 5* under the SU(5) GUT gauge group, the branching 
ratio of the radiative decay is given as B 1 ~ 5 x 10~ 6 . Even with such a small _B 7 , extra 
photons produced at a late time are cosmologically dangerous. 
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First, the saxion decay is constrained from the diffuse x-ray background spectrum. When 
the saxion has a lifetime longer than the cosmic time of the recombination {j„ > 10 12 sec), 
the produced photon directly contributes to the spectrum, and then a very stringent upper 
bound on x (p a /s) comes from the present observation of the spectrum. (The details are 
found, for example, in Ref. [T5|].) 

If the saxion lifetime is in between 10 6 sec and 10 12 sec, extra radiation energy produced 
by the saxion decay may alter the cosmic microwave background spectrum from the black- 
body one. The observation by the COBE satellite [[U| gives the following upper bounds on 
the saxion abundances as 



5 7 — < 2.5 x KT 5 T D for 10 6 sec < r a < 10 10 sec, 
s 



B~— < 2.3 x KT 5 T D for 10 10 sec < r CT < 10 12 sec, 
s 



where To is the cosmic temperature at the saxion decay and estimated as 

2/3 

M, 



T, 



D 



0.60 



G 



Pa 

S 



(23) 
(24) 



(25) 



These constraints on the saxion abundance are summarized in Fig.|| for the case B 1 = 
5 x 10~ 6 . Note that the constraints from the x-ray background spectrum and the CMBR 
become weaker for the suppressed I? 7 . We can now compare these cosmological constraints 
with the saxion abundance estimated in the previous section. 



B. Case 1: Saxion trapped at the origin 

First, we consider the case that the saxion is thermally trapped at the origin after 
the primordial inflation. For the case that the reheating temperature of the inflation is 
higher than the critical temperature (Tr > Tc ~ 100 GeV), the saxion abundance [Eq([L^)] 
significantly exceeds the upper bound Eq. ([2~T|) from the big bang nucleosynthesis. To avoid 
this the saxion should decay before it, however the lifetime of the saxion is longer than 1 
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sec if m a ^ 100 MeV. Even for the case m a ~ 100 MeV, the energy of the decay produced 
axion easily exceeds the bound Eq.([H|). Therefore, in this case, we conclude that Tr < Tq. 

The saxion abundance [Eq.([D^)] for Tr < Tc takes its minimum value at the lowest 
Tr = 10 MeV, since the dilution of the inflation becomes most effective. In Fig.|l] we show 
this lower limit of the saxion abundance with various cosmological constraints. Note that now 
for the case Tr = 10 MeV, the PQ scale as high as 10 16 GeV is allowed in the consideration 
of the relic abundance of the axion oscillation. It is found that all the region for the saxion 
mass is cosmologically allowed if we take the sufficiently low reheating temperature Tr ~ 10 
MeV. Furthermore the stable saxion with mass 10 eV ^ vn a < 10 keV will constitute a dark 
matter of our universe, since the saxion could achieve Q a = (p a /p c )o = 1 without conflicting 
with the cosmic x-ray background observation. 

In the above result an extremely low reheating temperature is crucial to dilute the saxion 



(and also the gravitino). The new inflation model |T7j gives such a situation. For example, 



Ref. |18| gives a new inflation model which solves the initial condition problem, where the 
reheating temperature is given by 

T R ~ m^MT, (26) 

where the gravitino mass should be larger than about 50 keV since Tr ^ 10 MeV. With this 
reheating temperature the saxion abundance fllSD is also shown in Fig.|l|. It can be seen that 



with this new inflation model |18j the saxion with m a ~ 50 keV-1 MeV survives the various 
cosmological constraints. 



C. Case 2: Saxion displaced far from the origin 

Next we turn to the case that the saxion is displaced far from the origin after the 
primordial inflation. In order to avoid the overproduction of the gravitino we assume the 
reheating temperature of the inflation to be lower than T a (even for the case that 1713/2 — 
m a < 1 keV). The saxion abundance for T R < T a is given as Eq.flTEl) and takes its minimum 
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value when the minimum Tr = 10 MeV. In Fig.[| we show the cosmological constraints on 
the initial amplitude <To in the saxion abundance Eq.(|T8|) with Tr = 10 MeV. Note that in 
this case the upper bound on the PQ scale is fpQ ^ 10 16 GeV. We find that most of the 
saxion mass (or the PQ scale) is cosmologically allowed if we take do ~ Jpq. And when the 
saxion has a heavier mass m a ~ 10-100 MeV, the initial amplitude can take a value from 
°o ~ fpQ to cr ~ Mq. Similar to the previous case, the saxion with mass 10 eV < m a < 10 
keV could be a dark matter of our universe. 

To summarize, the hadronic axion model we proposed is cosmologically viable in the 
both cases as far as the reheating temperature of the inflation is low enough. 

To end this section, we emphasize that it has been assumed that no extra entropy produc- 
tion takes place other than the primordial inflation. If such a late-time entropy production 
like a thermal inflation [H| or an oscillating inflation |l^,^,^l|] occurs, tremendous entropy 



is released at a late-time in the history of the universe (but at least before the big bang 
nucleosynthesis) and the saxion abundance is diluted by about 10 10 -10 20 and cosmological 
constraints on it are extensively relaxed. Thus our model becomes cosmologically viable 
even with the primordial inflation model with a high reheating temperature. 

IV. LATE DECAYING SAXION AND STRUCTURE OF THE UNIVERSE 

As discussed in the previous section, the cosmic energy of the saxion oscillation is abun- 
dant compared to the radiation and it may decay at very late time of the universe. In this 
section, we would like to argue that this feature of the saxion particle could explain well the 



present observed power spectrum of the cosmic structure with some cold dark matter p2 
i.e., the saxion can be a candidate for a late decaying massive particle (LDP) |23| , |24|| . 

It is well known that the predicted spectrum of the standard cold dark matter (CDM) 
with Qq — 1 and h ~ 0.7 could not fit well all the observational data simultaneously. The 
cosmic background explorer (COBE) observed the anisotropy of the temperature of the 
cosmic microwave background radiation (CMBR) [EH] and gave the normalization of the 
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spectrum at a large scale (A ~ 10 3 /i _1 Mpc). However, this normalization predicts too much 
fluctuation at a small scale (A ~ 10/i _1 Mpc) than galaxy distribution observations ||26| . The 
late decaying massive particle (LDP), once dominating the energy of the universe before its 
decay, can delay the epoch of the matter-radiation equality compared to the standard CDM 
model, which results in the preferred spectrum 

Ref. |24f argued that the saxion can be accounted to be the LDP, if the saxion energy 



density satisfies 



P - > n x El t (27) 
s s 



to cause the saxion domination era, and 

~ xM G , (28) 



2 



Pa 

S 



to explain the observed power spectrum. Here x is defined as 

* = o. a [(M) _xj . (29) 

In the following we will take x = 2.2, which corresponds to the case that Q = 1 and 
h = 0.7. Note that the saxion LDP also survives the big bang nucleosynthesis constraint 
(EHl). In Fig.||| we show the saxion abundance required to become the LDP by using its 



lifetime [Eg. (p0|)1. In the same figure, we also show various cosmological constraints. We 
find that the saxion whose mass is 

lOMeV > m a > lOOkeV (30) 

can be the LDP. However, such a saxion is severely constrained by the CMBR constraints 
[Eq. -P°| When the saxion becomes the LDP, the CMBR constraints give the upper 
bounds on the radiative branching ratio as 



Note that the saxion with mass ~ 10 MeV is free from the CMBR constraint since the lifetime 
is short. 
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5 7 < 8.8 x 10~ 6 , for 10 6 sec < r CT < 10 10 sec, (31) 

5 7 < 8.2 x 10~ 6 , for 10 10 sec < r CT < 10 12 sec. (32) 

Therefore B 1 = 5 x 10~ 6 for the simple case that the PQ quark multiplets Qp and Q P are 5 
and 5* of the SU(5) gauge group is marginally allowed. However, B 1 depends very much on 
their representation. Here we claim that the observational bound of the CMBR at present is 
comparable to the prediction of the late decaying particle. In other words the late decaying 
saxion scenario presented here will be tested by more precise experiments of the CMBR in 
near future (e.g. the experiments by the PLANCK and MAP satellites). 



A. Axion CDM and Saxion LDP Scenario 

We will explain that the saxion abundance required for the LDP [Eqs.(|27|) and (f28[)1 
can be naturally offered in the present model. First, we consider the case that the saxion is 
trapped at the origin due to the thermal effects after the primordial inflation. Using the relic 
abundance of the saxion Eq. flT5|) , the condition for the LDP [Eq.(p8|)1 leads to the required 
reheating temperature of the inflation asQ 

Tr - 23 11/10 te) w • 

K !30 MeV f^V 2/5 (J^Y 2 (^^f (^L_r, (34) 
\fp Q J viMevy VIOOGeVy v 10 GeV / 

with the saxion mass Eq. (|30|) . 

It is worth noting that if we take the critical temperature as Tq > 100 GeV, the reheating 

of the inflation, which induces the required saxion abundance of the LDP, completes before 



11 For example, a new inflation model by Ref. [18] with the reheating temperature (26) leads to 
the late decaying saxion which mass is 

m, * ,n 3/2 K 0.92 MeV fet) ' f (^)" , (33) 



which just lies in the required mass region [Eq. (p0|)1 . 
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the QCD phase transition. In this case the axion can constitute the dominant component 
of the dark matter of our universe if fpQ ~ 10 n -i0 12 GeV. This PQ scale corresponds to 
the saxion mass m a ~ 100 keV-1 MeV. Therefore, in the considering model, if one takes the 
gravitino mass 777.3/2 — m a ~ 100 keV-1 MeV, the axion can be the CDM of our universe 
and moreover the saxion can be the late decaying particle to cure the difficulty in the power 
spectrum of the structure of our universe in the axion CDM model. 

Even if the saxion is displaced from the origin after the inflation, a similar argument 
holds. Using the saxion abundance for T a > Tr [Eq. (|i~8|) l, the saxion which mass in Eq.([30"l) 
becomes the LDP when 

— ~ 0.91 1 „? , 



1-5 xlO 3 -^- — ^_ • (35) 



fpQ T R p 

m a \ 9 / 4 / T R y 1 ' 2 ( f 
lMevJ UOOGeVy yiO 4 GeV / 

Therefore if we take the reheating temperature as Tr ^ Aqcd, the gravitino mass 7773/2 — frier 
100 keV-lMeV also leads to the axion CDM. Thus the axion CDM and the saxion LDP 
scenario is achieved with the moderate Tr and a Q . 

B. Moduli CDM and Saxion LDP Scenario 

So far we have assumed that there exists no extra dilution mechanism other than the 
primordial inflation in the history of the universe. Here we will briefly discuss the case 
that there exists the late time entropy production such as the thermal inflation |19|] or 



the oscillating inflation |13| , pO| , pT| . The idea of the late time entropy production is that 



a mini inflation takes place at a late time of the universe (at least before the big-bang 
nucleosynthesis) and produces tremendous entropy at the reheating epoch, which dilutes the 
unwanted particles extensively. If one assumes this dilution mechanism, the cosmological 
constraints on the saxion abundance become much weaker even for the saxion abundance 
[Eq.flTTD] with <t ~ M G . Here we would like to stress that the saxion can play the role of 
the LDP even in this case. 

16 



In the universe with the late time entropy production it is worth noting that the string 
moduli could be the CDM.Q The moduli are scalar fields associated with flat directions which 
characterize classical vacua of the string theory and have only gravitationally suppressed 
interaction. The effects of the SUSY breaking give the moduli masses comparable to 
the gravitino mass (m^ ~ WI3/2) [ 30 1 - In fact if the moduli mass is < 200 keV, the 



condition f2^ ~ 1 could be achieved without conflict to the X(7)-ray observations after 



an appropriate amount of late time entropy production is taken into account 11151 , 1271 . 128 
Note that it extensively dilutes various dark matter candidates as well as the unwanted 
particles and they could not give a sizable contribution to the energy density of the universe. 
Therefore, the string moduli seem to be a well motivated candidate for the CDM in the 
universe with the late time entropy production.^ 

When we consider the moduli CDM, we find that the saxion can play a role of the LDP 
if the saxion is displaced at 07 ~ Mq after the primordial inflation.^ Note that the moduli 
(ft is also expected to be displaced from the true vacuum at the order of the Planck scale, 
i.e., ~ Mq- The both fields begin to oscillate at almost the same time, when the Hubble 
parameter becomes H ~ ~ m a (~ m 3/ / 2 ). Using the initial amplitude of the moduli 
(saxion) oscillation <f> (cr ), the energy densities of the oscillations are related as 

P, ^ H {^j • (36) 

This equation also holds after some late time entropy production, since the both energy 
densities are diluted at the same rate. Therefore, the saxion abundance is given by 



12 Cosmology of the string moduli in the gauge-mediated SUSY breaking is discussed by Ref. 
[13,27, 2^] in the universe with the thermal inflation, and by Ref. [13,2fj,2£] in the universe with 
the oscillating inflation. 

13 The cosmic axion with relatively high decay constant fpQ ~ 10 16 GeV may be another candidate 
for the dark matter 31,8|||. 



14 This initial condition is easily achieved by the effects of the supergravity. 
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This abundance naturally explain the conditions of the saxion LDP Eqs.(p7|) and (|28|). In 
fact, Eq.fl2"E|) leads to the ratio of the initial amplitudes 

a _ 0.32 mfMj /4 
To^n^f(p c /s ) 1/2 

5-4 / m a \ s/ 4 ( f \ 1 ,„ , 

-Q^KlMeVj ll0 4 Gevi ' 1 } 



with the saxion mass in the region fl3"D|). Here the moduli CDM implies ~ 1. It should 
be noted that both initial amplitudes is expected to be of the order of the Planck scale, 
cto ~ 0o ~ Ma- Thus for Vt^ ~ 1 the saxion with the initial amplitude gq which is slightly 
larger than <f) ~ Mq could be the LDP even with the late time entropy production in the 
history of the universe, if the saxion has a mass in the region (|30|) . On the other hand the 
string moduli could be the CDM of our universe (fl^ ~ 1) if its mass is less than about 100 
keV. Therefore, if m a ~ ~ 77^3/2 ~ 100 keV, the string moduli could be the CDM of our 
universe while the saxion could plays a role of the LDP so that the present observed power 
spectrum is well explained. 

It is interesting that this scenario will be tested in various cosmological experiments. The 
moduli CDM with mass ~ 100 keV will be tested in the future by the x-ray background 
observations with high energy resolutions The line x-ray spectrum from the moduli 

CDM trapped in the halo of our galaxy will be seen. And the future CMBR observations, 
as mentioned before, reveal the saxion LDP possibility. 



V. CONCLUSIONS 

In this paper, we closely investigated the cosmology of the hadronic axion model in the 
gauge mediated supersymmetry breaking, which we had proposed in a previous paper . In 
particular we focused on the implications given by the saxion which is the scalar component 
of the axion supermultiplet. The saxion has a tiny mass comparable to the gravitino mass 
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and thus its life time tends to be long and decay after the primordial nucleosynthesis. 
Furthermore since the saxion is a scalar field, it will be abundant in the early universe in 
the form of the coherent oscillation. Thus the saxion may be cosmologically harmful. 

We estimated the energy density of the saxion in two distinctive cases for the cosmological 
evolution of the saxion coherent mode. Then we examined whether or not the saxion survives 
the various cosmological constraints from the nucleosynthesis, from the critical density limit, 
from the diffuse X-ray background spectrum, and from the cosmic microwave background 
spectrum. We showed in both cases that the model we proposed is cosmologically viable, 
namely it satisfies the severe cosmological constraints on the saxion field, if the reheating 
temperature of the primordial inflation is sufficiently low. 

Furthermore we argued that the saxion can play a role of the late decaying particle in 
the structure formation of the universe, giving a better fit of the power spectrum of the 
density perturbation than the standard cold dark matter model. 
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FIG. 1. The lower bound of the saxion abundance [Eq,(|i~5|)1 for the case Tq = 100 GeV and 
/ = 10 4 GeV (the thick solid line). We take the value of the radiative branching ratio -B 7 = 5 x 10 -6 . 
The upper bounds from the various cosmological constraints are also shown. For the saxion mass 



larger than about 50 keV, the predicted saxion abundance in the new inflation model [18] is shown 
with the thin solid line. 
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FIG. 2. Various cosmological upper bounds on the initial amplitude of the saxion oscillation 
[Eq.(|lD] for the case T R = 10 MeV (< T a ). We take / = 10 4 GeV and B 7 = 5x 1(T 6 . We also 
show the initial amplitude ctq = fpQ by the solid line. 
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FIG. 3. The abundance of the saxion which is required to become the LDP (the solid line). Here 
we take / = 10 4 GeV, S7o = 1 and h = 0.7. The upper bound from the big bang nucleosynthesis 
and the lower bound from the saxion domination are denoted by the dotted lines with shadows. 
The dot-dashed lines represent the upper bounds of the saxion abundance from the CMBR for the 
cases that the radiative branching ratio B 7 = 5x 10~ 6 . 
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